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abstract 

The structure of moving nonlinear excitations in one-dimensional electron- 
phonon systems is studied semi-phenomenologically by using an effective ac- 
tion in which the width of the nonlinear excitation is treated as a dynamical 
variable. The effective action can be derived from Su, Schrieffer and Heeger's 
model or its continuum version proposed by Takayama, Lin-Liu and Maki 
with an assumption that the nonlinear excitation moves uniformly without 
any deformation except the change of its width. The form of the action is 
essentially the same as that discussed by Bishop and coworkers in studying 
the dynamics of the soliton in polyacetylene, though some details are differ- 
ent. For the moving excitation with a velocity v, the width is determined by 
minimizing the effective action. A requirement that there must be a mini- 
mum in the action as a function of its width provides a maximum velocity. 
The velocity dependence of the width and energy can be determined. The 
motions of a soliton in polyacetylene and an acoustic polaron in polydiacety- 
lene are studied within this formulation. The obtained results are in good 
agreement with those of numerical simulations. 
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1 Introduction 



It is generally accepted that dynamical behaviors of nonlinear excitations in 
one-dimensional electron-phonon systems such as conjugated polymers play 
an important role in determining physical properties of those systems. [0 
The motions of the nonlinear excitations have been studied by analytic meth- 
ods 0, U and by numerical simulations. [|, |], H I n recent papers, || ||, [10 



TT| , |l"2|| using Su, Schrieffer and Heeger's (SSH) model |H| modified to include 



an external electric field, we have investigated the dynamics of a charged soli- 
ton in polyacetylene by numerical simulations in the presence of an electric 
field. We have found the following results: i) The soliton has a saturation 
velocity of about four times the sound velocity. || |[ ii) The width of the 
soliton decreases as the soliton velocity increases. || Furthermore, the time 
dependence of the width determined from the spatial variation of the bond 
order parameter involves an oscillatory component, which is related to the 
amplitude mode around the soliton. [0 iii) The overall behavior of the total 
energy of the system as a function of the soliton velocity v is well expressed 
as 

Ae tot (t;) = -(M s ^/2)ln(l-^/<) (1) 

with As tot (v) = e to t(v) — e to t(0), where v m is the maximum velocity and M s 
the soliton effective mass of about four times the bare electron mass. |TIJ 

Bishop and coworkers |J also have studied the dynamics of the soliton 
with numerical simulations and shown the saturation of the soliton velocity 
whose value is in rough agreement with our numerical result. Furthermore 
they have given a phenomenological explanation for the saturation mecha- 
nism that the time necessary for an ion group on a lattice point to respond 
to the reverse of the dimerization pattern caused by the soliton translation 
can not be shorter than the characteristic time of the lattice vibrations, ujq . 
Here ujq is the renormalized optical phonon frequency. Alternatively they 
have explained the saturation by another independent method. According 
to that method, the width £ of the soliton moving with a velocity v is de- 
termined by minimizing an effective action. A requirement that the action 
must have an minimum as a function of £ provides the saturation (maximum) 
velocity. 

In this paper, we reconsider the phenomenological method for describing 
the structure of moving nonlinear excitations in one-dimensional electron- 
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phonon systems by using the effective action originally proposed by Bishop 
and coworkers. This method is applicable not only to the soliton in poly- 
acetylene but also to a polaron or other nonlinear excitations in conjugated 
polymers. The effective action is made up of a potential energy part and a 
kinetic one, and the width of the nonlinear excitation is regarded as a dy- 
namical variable. The potential energy is defined by the energy of a static 
nonlinear excitation for each value of the width. Both parts are calculated 
numerically or analytically by using the SSH model [^] or its continuum 
version, e.g. Takayama, Lin-Liu, and Maki's (TLM) model in the case of 
the soliton in polyacetylene. [II] Using the effective action, we determine the 



saturation velocity, the relation between the velocity and the width, and the 
velocity dependence of the energy of the nonlinear excitations. We concen- 
trate on the dynamics of a soliton in polyacetylene and of an acoustic polaron 
in polydiacetylene. |T5| The latter is formed by interactions between a free 
electron and acoustic phonons, and its dynamical behavior has been investi- 
gated analytically by Wilson. || The width and the energy of the acoustic 
polaron have been presented as the function of its velocity. For both cases, 
we show that the dynamical behaviors evaluated from the action are in good 
agreement with those obtained from the numerical simulations. Furthermore, 
for the case of the acoustic polaron, it is indicated that the results derived 
from the effective action are consistent with the analytic results by Wilson. 

This paper is organized as follows. In § 2 we define the effective action for 
the moving soliton and its structure is analyzed. The results are compared 
with those obtained by the numerical simulations. The same analyses are 
performed for the acoustic polaron in § 3. We present conclusions and discuss 
the results in the final section. 



Structure of Moving Soliton in Peierls-Dimerized 
Linear Chain 
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2.1 Effective action 



In order to discuss the effective action describing the behavior of a moving 
soliton, we start with the SSH model, 

#ssh = ~X)(*o - uyn)[c ] na c n+l(7 + h.c] + —Y^vl + -o-IX' ( 2 ) 

no- ^ n 

where y n = w n+ i — u n . The quantity £ is the nearest neighbor transfer 
integral of 7r-electrons in the regular lattice, a the electron-lattice coupling 
constant, and u n the lattice displacement of the n-th site from its equidistant 
position. The operators cj^ and c na create and annihilate a 7r-electron with 
spin a at the n-th site, respectively. K is the force constant mainly due to 
the cr-bond, M the mass of a CH unit. 

We consider a uniformly moving soliton with a velocity v. The soliton 
is regarded as a spatially extended object with an internal degree of free- 
dom characterized by its width. The lattice displacement at the n-th site is 
assumed to be expressed as, 

u n = (-) n u tanh( j, (3) 

where u denotes the amplitude of the staggered lattice displacement in the 
perfectly dimerized state, a the lattice constant, and £ the soliton width. 
Here, in a collective-coordinate theory, the center of mass of the soliton and 
its width are treated as dynamical variables characterizing the structure of 
the moving soliton. |TB, [T7], P| These two collective coordinates depend on 



time. In the present formulation the center of the soliton moves with the 
velocity v and the width is assumed to be determined by the velocity. We 
ignore any other internal excitations such as the amplitude mode. 

The effective action is defined by the potential energy minus the kinetic 
energy, 

A(r;v) = V(r)-—, (4) 

with 

r = M, ^Mto 2 , ( 5) 

£o 3£ \v s J 

where £ is the width of a stable static soliton and v s the sound velocity 
(= a^J K/M). V(r) denotes a potential energy which is given by the sum of 
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the electronic energy and the lattice elastic energy of the system containing 
a soliton with the relative width r. The second term on the r.h.s. of eq. (^) 
comes from the lattice kinetic energy The details of the derivation of the 
effective action are given in Appendix. In Fig. 1 the potential energy V(r) 
is shown, where the optimized energy at r = 1 is taken to be zero. We can 
neglect the constant part independent of r because in the present method it 
is essential to find a minimum of the action as a function of r. 

The action can also be defined from the continuum version of the SSH 
model, i.e. the TLM model. In this case V(r) is essentially identical with 
the creation energy of a static soliton for the fixed relative width r, and the 
coefficient of the kinetic energy term is expressed as 



HV o \ 2 2 ' \ / 



where 2A is the Peierls gap, A(= a 2 /Kt ) the dimensionless electron-phonon 



coupling constant, vp the Fermi velocity and uJq(= yAK/M) the bare optical 
phonon frequency. It is straightforward to confirm that f3 v in eq. (||]) reduces 
to that in eq. (§) in the continuum limit ( or in other words, the weak 
coupling limit). 

Here we note that, in the study by Bishop and coworkers, || V(r) is 
approximated by a quadratic function of r which coincides to V(r) calculated 
from the TLM model at two points r = and 1. Its explicit r-dependence 
is shown in Fig. 1 by a dashed curve. Then the effective action is rewritten 
as, ]T9| 



A=-A + A (l--)(l-r) 2 -^, (7) 

7T 71 r 

The first term on the r.h.s. is the energy of the soliton at rest. However this 
approximation is too crude for our purposes as will be seen in the following. 



2.2 Saturation of velocity 

The maximum velocity is determined from the requirement that the action 
must have at least one minimum as a function of r. If the approximation eq. 
(|?D is used, the requirement is fulfilled when (3 V < /3 C = (2/3) 3 (l — (2/tt)). 
For f3 v = /3 C , the minimum point of the action is located at r = 2/3. From 
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(3 C} the maximum velocity is estimated to be 



2/2 

v m3X = ~\ir\(l )v F u Q /A ~ 3.74v s (8) 

O V 7T 

This velocity is appreciably larger than the saturation velocity (~ 3.22v s ) 



obtained from the numerical simulation in the previous works. [10 

Next we use the exact expression of V(r) defined from the SSH model in 
order to calculate /3 C , the critical value of (3 V . (3 C is defined as a value where 
the minimum of the derivative of the action with respect to r becomes zero, 
and is calculated numerically. In Fig. 2 typical examples of the action and 
its derivative are shown for (3 V < (3 C (a), (3 V = (3 C (b) and (3 V > (3 C (c). From 
P c , the maximum velocity is evaluated to be t> max = 3.18t> s , which is in good 
agreement with that obtained in the numerical simulation. When the action 
defined from the TLM model is used we obtain also a rather good result, 

"max 3.22f s . 

2.3 Velocity dependence of the soliton width 

The width of the moving soliton is determined from the condition of min- 
imizing the effective action, as described above. The velocity dependence 
of the width obtained by using the effective action defined from the SSH 
model is shown in Fig. 3 along with the result of the numerical simulation 
for comparison. For the numerical result we used the width calculated from 
the bond order parameter. || 

The relation between the width and the velocity obtained from the action 
is in good agreement with that of the numerical simulation. We fit the 
relations to the following form through the least square fitting in the region 
v < 2v s , 

and obtain for the analytic curve, 

a ana i y = 2.4 x KT 2 , 

&analy = 5.3 X KT 4 , 

and for the curve of the numerical simulation, 

Onum = 1-6 X 1(T 2 , 



bnum = 9.7 XlO" 4 
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The small difference between the analytic and numerical results may be 
due to the inaccuracy in the velocity and the width calculated from the lat- 
tice displacement pattern in the numerical simulation. These quantities have 
a rather poor accuracy due to the tail structures in the lattice displacement 
pattern after the soliton has traveled a rather long distance. The tail struc- 
tures may particularly give an inaccuracy to the width. Taking account of 
this inaccuracy, the agreement between the semi-phenomenological analysis 
and the numerical simulation seems to be satisfactory. 



2.4 Velocity dependence of energy 

The energy of a moving soliton is represented by 

E{v) = V{r) + Sl, (10) 

where the relative width r is the function of the velocity v as shown in 
Fig. 3. Figure 4 indicates the velocity dependence of the soliton energy, in 
which results of the numerical simulation are also shown for the total and the 
lattice kinetic energies of the system containing the moving soliton. Although 
a slight difference in magnitude between the analytic result for the total 
energy and that of the simulation is seen, the functional form of the analytic 
result reproduces well the numerical one. Here the energies calculated in the 
simulation involve the contributions of the lattice vibrations (tail structure) 
emitted from the moving soliton; this contribution is not directly related 
to the energy of the soliton. If this effect is taken into account, the slight 
difference would become smaller. 

In our previous paper, we have shown that the energy of the moving 
soliton looks to diverge logarithmically as a function of the velocity (eq. ([!])). 
Here we deal with the behavior of the energy around the maximum velocity, 
using eq. (|i~0f) . The derivative of E(v) with respect to v is expressed as 

dE(v) dr (dV (3 2 \ 2(3 
H v, 



dv dv \ dr r 2 ) r 

2(3 2 dr 2(3 

= 2 V T + — 

r z dv r 



2 V - + ^, (11) 



where (3 is a constant defined by (3 — (3 v /v . In deriving the last expression 
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we have used the relation between r and v, 



dA 
dr 



dV 
dr 



(3 



0. 



(12) 



The velocity derivative of r is obtained by differentiating eq. (112 ) by v, 



d_ 

dv 



v or J 



dr 
dv 



'd 2 V _ 2/3 
v dr 2 

iax 



-V 



r" 



0. 



(13) 



On the other hand, when v = v r 
that - = — and -^-r- satisfy following relations 



/3 C ), it is obvious from Fig. 2(b) 



dr 



dr 2 



dA 
dr 

d 2 A 



dr 2 



0. 

d 2 V{r r 



2(3 



dr 2 



0. 



(14) 
(15) 



From eqs. (13) and fll5]) the derivative of r with respect to v diverges to minus 

dE 

infinity at v = v m ax- Then — — diverges as v approaches v max . Noting that 

dv 

the value of r itself does not diverge at v — v mSLX , we can safely assume that 
r(v) = r(t> max ) + C (v — f max ) M with < \i < 1 near f max . This assumption 
will be confirmed by the result. If we put r(v) = r(t> max ) + Ar, then Ar can 
be estimated from eq. (fi2|) in the following way. 

dA 

We expand the expression for — — in terms of Ar and Av(= v — v max ) up 

dr 

to second order. 

1 d 3 A 



dA 
dr 



fmin(f) 



2 dr 3 

+4- 

0, 



fminf^n 

2 



(Av 



(Ar) 2 + 





2^ D 



mm V "max 



)] 



-Av 



-(Ar)(Av) 



(16) 

where we have used eqs. (14) and ([15]). It is obvious from what we have 
assumed for Ar that the last two terms in the above expansion are irrelevant. 
It is also clear from the general behavior of A(r) that 

» A > 0. (17) 

^*min(^max) 



dr 3 
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As a result, we end up with 

Ar OC V^max - v, (18) 

consistently with the first assumption. Therefore we can conclude that 

E(v) ~ E(v max ) - Cy/v max - v, (19) 

by combining the above result for Ar with eq. ([0]). E(v max ) itself is finite 
dE 

but — — diverges at v — v max . This singularity would have been seen as a 
dv 

logarithmic singularity when we analyzed the overall energy- velocity relation 
obtained from the numerical simulation. 



3 Dynamics of Acoustic Polaron 
3.1 Effective action 

An acoustic polaron can be described in terms of the SSH Hamiltonian 
by postulating that the system has only one electron. In Fig. 5 the self- 
consistent solution is shown for the charge distribution, the lattice displace- 
ment and the bond order parameter, respectively. The values of the pa- 
rameters used are t = 0.64eV, a = 0.38eV/A, K = 2.4eV/A 2 , a = 4.9A, 
uq = 0.96 x 10~ 2 eV, which might be appropriate for polydiacetylene. The pe- 
riodic boundary condition is assumed. In an infinite system, a static acoustic 
polaron can be well expressed by the form, || 

Un f na\ , . 

u n = -ytanh(— j, (20) 

h 2 Mv 2 

with £ (= t4t) the width of a stable acoustic polaron and u (= 

2aa 2 / Mv 2 ) the amplitude of the lattice displacement at a long distance from 
the center in the case of a stable polaron, || where we assume the center to 
be at the origin. For a stable polaron the relative width r is equal to unit. 
The prefactor in eq. (p0|) is determined so as to minimize the total energy 
of the system when the polaron width £ r is fixed. The change of the lattice 
spacing affects directly the transfer integral of the electron, 

auo 2 / na \ , N 

y n = u n+1 -u n ~ -— ^sech I — I . (21) 
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Figure 6 indicates the acoustic polaron in the infinite system expressed by u n 
and y n . In the numerical simulations of the motion of the acoustic polaron, 
we use a finite size system and impose the periodic boundary condition. This 
condition leads to additional terms in eqs. ( p0|) and (|21|). 

Un = _^tanh(^Uc™> (22) 
r \£ rJ 

fc = -^d'fjlUc, (23) 



with 

au ^ t.2^ a \ 2u o / 0/A 

G= w^J N/ r h krr)^N;- < 24 > 

A uniformly moving acoustic polaron will be described in the u n repre- 
sentation as 

Uq (na — vt\ C , . /ri _. 

u n = tanh — H (na — vtj. (25) 

r \ £ r J a 

Inserting this expression into the kinetic energy part in the SSH model, we 
get 

M_. 9 M „ „ f i -/z-uA 2 ) 2 



-Ek = — > w„ - — ^ «n / < — -sedr ^ 

2 V 2a °7-iVa/2 \e r 2 V & J Nar I 



2Mv 2 u 2 ( 1 1 



3^o^ 3 Nar 2 



(26) 



The effective action for the moving acoustic polaron is defined similarly to 
the case of the soliton as, 



The potential term V(r) is calculated by substituting eq. (|23|) into the SSH 
model in the same way as for the soliton in polyacetylene. We show the 
potential V(r) in Fig. 7. 
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3.2 Saturation of velocity 

The saturation of the velocity in the case of the acoustic polaron has been 
discussed by Wilson and estimated to be nearly equal to the sound veloc- 
ity. || The dynamical simulation of the acoustic polaron in the presence of 
an electric field shows that it is about 0.8 ~ 0.9v s . pOf By using the effective 



action it is found that the acoustic polaron becomes zero width at the sound 
velocity as is shown in the next section. Then the maximum velocity is ex- 
pected to be the sound velocity from this fact. The reason why it is slightly 
smaller in the numerical simulation will be that the computing time is finite 
and also that the lattice vibrations excited by the motion of the polaron may 
be acting as scatterers yielding a small resistance against the motion. 

3.3 Velocity dependence of width 

The solid line in Fig. 8 indicates the relation between the width and the 
velocity of the acoustic polaron determined from the action. The data of the 
numerical simulation is also shown, pi The agreement between the analytic 



result [gl] an d the numerical one is quite good. It seems that the width of 
the acoustic polaron becomes zero at v = v 8 . According to the analytical 
study by Wilson, [|| who used instead of the effective action the equation 
of motion for the lattice displacement and the Schrodinger equation for the 
electronic wave function in the continuum limit, the velocity dependence of 
the width is of the form 

(28) 

which is indicated by the broken line in Fig. 8. Furthermore it has been 
shown that the electronic energy and the lattice potential energy as well as 
the lattice kinetic energy are expressed as the functions of r as|| 

(29) 
(30) 

-) 2 - (3D 





U 


e e = 






2U 


CLP = 


3 r 3 




2U 


^LK = 


3 r 3 
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Here 



_ m 8«V 



If we use these expressions, we obtain the same result as eq. ( p8"D with the 
help of the effective action, which supports the appropriateness of the semi- 
phenomenological method. In fact it is straight forward to confirm that the 
lattice kinetic energy eq.(^) reduces to eq. fl3ll) in the infinite system. 

Here note that we recognize an appreciable difference near the sound 
velocity between the velocity dependences of the widths calculated from the 
action and from the continuum approximation by Wilson in Fig. 8. This 
may be due to the discreteness of the SSH model, since the width of the 
acoustic polaron becomes quite narrow near the sound velocity. 

3.4 Velocity dependence of energy 

The energy of the moving acoustic polaron is given by 

E(v) = V(r) + E K . (33) 

The velocity dependence of the energy E{v) and that of the kinetic energy 
part -Ek are shown in Fig. 9 together with the results of the numerical 
simulations. Both of the energies diverge as the velocity approaches the 
sound velocity. The agreements between the analytic and numerical results 
are excellent. These results agree also with those obtained by Wilson. This 
fact along with the results in the previous subsections indicates that the 
present method in terms of the effective action is essentially equivalent to 
Wilson's analytic method and that both methods can reproduce the results 
of the numerical simulation quite well. 



4 Summary and Discussion 

We studied the structure of a moving soliton and a moving acoustic po- 
laron in the one-dimensional electron-phonon system by means of a semi- 
phenomenological method by introducing the effective action. The effective 
action for a moving nonlinear excitation consists of the potential energy with 
its width as the "coordinate" and of the kinetic energy. In the continuum 
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limit (the TLM model) for the soliton problem, the potential energy is ex- 
pressed by the creation energy of the soliton with a given width. It is assumed 
that the width can be regarded as a dynamical variable which depends on 
the velocity. The widths of the moving nonlinear excitations are determined 
so as to minimize the effective action. The requirement that the action must 
have a minimum gives the upper limit of the velocity, and its value is found 
to be nearly the same as the saturation velocity obtained by the numerical 
simulation. The velocity dependence of the width and that of the energy 
calculated by the semi-phenomenological method present almost the same 
behaviors as the results obtained from the dynamical simulations. We may 
reasonably conclude that the dynamical behaviors of the soliton and the 
acoustic polaron can be described by the phenomenological effective action. 

From the analysis using the effective action, it is found that for the case 
of the soliton there occurs the divergence not in the energy itself but in its 
derivative with respect to the velocity at v — v max . In contrast to the case of 
the soliton, the energy of the acoustic polaron diverges at v — v s . 

The dynamical behavior of the acoustic polaron is described quite well 
by the effective action while such a description was not much accurate in the 
case of the soliton in polyacetylene. One of the reasons for this is the acoustic 
polaron has no extra internal excitation in contrast to the case of the soliton 
which has several localized modes omitted in the present treatment. 

In this paper we have focused our attention on the soliton and the acoustic 
polaron. The semi-phenomenological method developed in this work is appli- 
cable to the cases of other nonlinear excitations in one-dimensional electron- 
phonon systems such as an optical polaron. We are now studying this subject, 
the result of which will be presented elsewhere. 
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We show the derivation of the effective action for the soliton. Although we 
treated the soliton width as a function of the velocity in the previous section, 
we consider it time depending variable here. 
Using eq. (|3|), the lattice kinetic energy is 

E M^ < = m\«l + t^e). (34) 



2 

The Lagrangian for the fixed velocity v is expressed as 

% a = **£W^_ m (35, 

with 



where V(£) denotes the potential energy which is given by the sum of the 
electronic and lattice elastic energies. Then the equation of motion is written 

(37, 

There can be the amplitude mode around the soliton, which leads to an os- 
cillation of the soliton width. Therefore, the equation of motion must have 
an oscillatory solution. If the equation of motion eq. fl3~7|) is solved, the am- 
plitude mode frequency may be evaluated as a function of the velocity, which 
we have studied numerically in the previous paper. We have confirmed 
that eq. ([37]) yields a reasonable value for the frequency of the amplitude 
mode in the case of the static soliton (v — 0). |22|] In the case of the moving 



soliton (v 7^ 0), however, the amplitude mode frequency derived from eq. 



(37) decreases with the velocity in contrast to the result of the numerical 
simulation where the frequency was an increasing function of v. The reason 
of this discrepancy is not clear at the moment. Presumably the above model 
will be too much simplified to discuss such detailed behaviors as delicate 
changes in inner degrees of freedom of the soliton. 

We have found in the numerical simulation that the amplitude mode is 



excited with the acceleration of the soliton. [F3] When the soliton moves 
with a constant velocity, the soliton width oscillates with an almost constant 
frequency and amplitude. In reality, however, there will be the damping 
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of the amplitude mode due to the interactions with other modes which are 
not included in the Lagrangian eq. fl35|) , and the amplitude of the width 
oscillation will damp and then the width settles down to a certain value, 
which depends on the velocity. The effect of the damping could be introduced 
by adding a phenomenological friction term in the equation of motion. After 
the elapse of sufficient time, the time derivative of the width and its second 
derivative can be both treated as zero. The equation of motion can be reduced 
to 

M 4fj - V(0) = 0, (38) 

which is the equation to determine the width of the soliton moving with the 
velocity v. The soliton width is determined by an extreme of the function 
inside the parentheses which we call the effective action. The effective action 
eq. (|4]) can be obtained by inverting the sign and introducing the relative 
width r in the function. 

The effective action for the acoustic polaron can be derived by a similar 
way as above. 



References 

[1] H. G. Kiess (Ed.): Conjugated Conducting Polymers, Springer Series in 
Solid-State Sciences 102 ( Springer- Verlag, Berlin, 1992). 

[2] G. M. Silva and Y. Wada: Synthetic Metals 43 (1991) 3713. 

[3] E. G. Wilson: J. Phys. C: Solid State Phys. 16 (1983) 6739. 

[4] W. P. Su and J. R. Schrieffer: Proc. Natl. Acad. Sci. 77 (1980) 5626; 
Phys. Rev. Lett. 46 (1981) 738. 

[5] E. J. Mele: Solid State Commun. 44 (1982) 827; Phys. Rev. B26 (1982) 
6901. 

[6] A. R. Bishop, D. K. Campbell, P. S. Lomdahl, B. Horovitz and S. R. 
Phillpot: Phys. Rev. Lett. 52 (1984) 671; Synthetic Metals 9 (1984) 
223. 



15 



[7] A. Terai: in Relaxation in Polymers, edited by T. Kobayashi (World 
Scientific Publishing Co., 1993), p. 269. 

[8] Y. Ono and A. Terai: J. Phys. Soc. Jpn. 59 (1990) 2893. 

[9] M. Kuwabara, Y. Ono and A. Terai: J. Phys. Soc. Jpn. 60 (1991) 1286. 

[10] Y. Ono, M. Kuwabara and A. Terai: J. Phys. Soc. Jpn. 60 (1991) 3120. 

[11] M. Kuwabara, Y. Ono and A. Terai: J. Phys. Soc. Jpn. 61 (1992) 2412. 

[12] M. Kuwabara and Y. Ono: J. Phys. Soc. Jpn. 62 (1993) 990. 

[13] W. P. Su, J. R. Schrieffer and A. J. Heeger: Phys. Rev. Lett. 42 (1979) 
1689; Phys. Rev. B22 (1980) 2099. 

[14] H. Takayama, Y. R. Lin-Liu and K. Maki: Phys. Rev. B21 (1980) 2388. 

[15] N. A. Cade and B. Movaghar: J. Phys. C: Solid State Phys., 16 (1983) 
539. 

[16] M. J. Rice and E. J. Mele: Solid State Commun. 35 (1980) 487. 
[17] M. J. Rice: Phys. Rev. B28 (1983) 3587. 

[18] M. J. Rice, S. Jeyadev and S. R. Phillpot: PR B34 5583 (1986). 

[19] This expression differs from what is obtained by expanding V(r) around 
r = 1 up to the second order of (r — 1). 

[20] Y. Arikabe, M. Kuwabara and Y. Ono: to be submitted. 

[21] Since the potential energy in the effective action is calculated numeri- 
cally in the case of the acoustic polaron, it may not be completely ap- 
propriate to call the result evaluated from the effective action "analytic 
result". However we use terminology "analytic" in order to emphasize 
that the numerical simulation is not necessary to analyze the structure 
of the moving nonlinear excitations. 

[22] A. Terai and Y. Ono: J. Phys. Soc. Jpn. 55 (1986) 213. 



16 



Figure Captions 

Fig. 1 The potential energy V(r) for the case of the soliton calculated from 
the SSH model (solid line). V(r) is equivalent to the creation energy 
of a static soliton as a function of its width. The energy is scaled by 
A (half of the Peierls gap) and the width by its ground state value £ - 
The deviation in energy from the ground state energy is plotted. The 
broken line indicates the approximation employed by Bishop et. al. || 

Fig. 2 The action (upper part) and its derivative by r (lower part) for 
(3 V < Pc (a), Pv = Pc (b), and (3 V > (3 C (c). In each case the potential 
energy (broken line) and the kinetic energy (dash-dotted line) are shown 
in upper part. 

Fig. 3 The velocity dependence of the soliton width determined from the 
effective action (solid line) and from the numerical simulations (broken 
line). The width £ is scaled by £ an d the velocity by the sound velocity. 

Fig. 4 The velocity dependence of the soliton energy determined from the 
effective action (solid line) and from the numerical simulation (broken 
line) . 

Fig. 5 Structure of the acoustic polaron for the case of the finite system 
size with periodic boundary condition. The charge density (top), the 
lattice displacement (middle) and the bond configuration (bottom) are 
shown as functions of the site number. 

Fig. 6 Structure of the acoustic polaron for the case of the infinite sys- 
tem size. The lattice displacement (top) and the bond configuration 
(bottom) are shown. The ordinates are in arbitrary unit. 

Fig. 7 The potential energy for the case of the acoustic polaron. 

Fig. 8 The velocity dependence of the width of the acoustic polaron deter- 
mined from the effective action £ A (solid line) and from the numerical 
simulation of the polaron motion £i at (solid dots). The broken line 
indicates Wilson's result £ w . 

Fig. 9 The velocity dependence of the energy of the acoustic polaron (solid 
line) and its kinetic energy part (broken line) determined from the 
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effective action. Those determined by the numerical simulation are 
indicated by dots and crosses, respectively. 
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